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Abstract 

The supergenerator of superparamagnetic system quantum evolution is in- 
vestigated in discrete orientation model (DOM). It is shown that the generator 
is J-self- adjoint one at the case of potential drift field agreed upon magnetic 
anisotropy of the sample investigated. Perturbation theory is used for spectral 
analysis. The qualitative dependence of resonance absorption spectrum on the 
relation between quantum and stochastic parameters is demonstrated. 
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Introduction 

Qualitative and quantitative description of relaxation Mossbauer spectra is based on 
choosing of some relaxation process model. In the most general model of motion of 
magnetic moment M(t) suggested by Brown [Ij] this motion is considered like diffusion 
process on radius \M\ sphere in a drift field stipulated by the magnetic anisotropy 
of the pattern. Under Born's approximation Mossbauer spectra line shape one can 
express in terms of resolvent of the generator of quantum evolution operator of a 
nucleus, averaged along all diffusion process M(t) trajectories. But computation 
of averaged evolution operator is connected with solving of complicated differential 
equation system of partial derivatives and can not be done in analytical form. To 
avoid this difficulty phenomenological models are used, where magnetic moment M(t) 
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position becomes discrete, and process of M(t) motion itself is replaced by Markovian 
process with final number of states, which is described by stochastic matrix with phe- 
nomeno logical transition probabilities of the magnetic moment M(t) from one state 
(of easy magnetization) to another (so called discrete orientation model (DOM)). In 
this case the problem of computation of averaged evolution operator obtains algebraic 
character and reduces to inversion of special matrices, which properties and order are 
determined by the number of easy magnetization directions, the structure of nucleus 
magnet electron system and by the symmetry of pattern. In this direction certain 
progress was reached for SP particles with cubic symmetry ||. Direct account of 
symmetry in frames of phenomenological DOM permitted authors M to reduce the 
calculation of absorption line shape of SP particles to inversion some matrices of 
8- and 6- order for cubic symmetry systems with negative and positive magnetic 
anisotropy constants respectively. 

For the case of large values of reduced barrier a = KVj kT , where K is magnetic 
anisotropy constant, k - Boltzmann constant, V - volume of SP particle, T - temper- 
ature, DOM parameters one can calculate by the lowest eigenvalues of Fokker-Plank 
equation ||, ||, giving evolution of M(t) in diffusion model 0, and line shape cal- 
culated in such a way is the leading term under a —>■ oo of line shape calculated on 
the base of diffusion model 0. For accurate calculation of the line shape error under 
substitution of diffusion model to discrete one (DOM) first of all detailed analysis 
of the generator of averaged evolution operator of a nucleus in frames of DOM is 
necessary. 

In present paper, developing symmetry considerations stated at 0, we investi- 
gate the structure of complete averaged evolution operator corresponding to DOM 
for cubic symmetry (with 6 and 8 easy magnetization directions on Fe 57 nucleus) and 
some others, maintain it's spectral analysis in frames of perturbation theory under 
assumption of slow relaxation and also describe qualitative picture of spectrum be- 
havior in general case (subject to temperature). In proper basis the generator matrix 
is composed of equal blocks of 8 x 8 older differing only by its diagonal elements 
composed by combinations of eigenvalues of the stochastic matrix giving DOM. Non- 
diagonal part of such block is determined only by the nuclear system structure and 
does not depend neither on stochastic variables nor on crystal lattice symmetry type 
(on magnetic anisotropy constant sign for cubic symmetry). 
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1 Basic Formulas 



In DOM magnetic moment evolution is denned by matrix P = S — I, where S is 
stochastic matrix of transition probabilities between easy magnetization states in unit 
of time. Vector g(t) of magnetic moment M(t) distribution is governed by ordinary 
differential equation 

~jr = p 9> 9 = (gi,92,---,9Ny, (i) 

where N - the number of process states (easy magnetization directions). In the case of 
cubic symmetry N is equal either 6 or 8 for positive or negative magnetic anisotropy 
constant respectively. Matrix P is naturally to call as Markovian one. 
Line shape is determined by expression 

oo 

tp(u) = -Re J e (i "- r/2)i Sp(V*G(t)Vp)dt, (2) 
o 

where u is the electromagnetic frequency, T - natural half-width, p - the initial state 
density matrix, V - the operator of interaction of the nucleus with electromagnetic 
field, G(t) - the result of averaging of the quantum evolution superoperator of the 
nucleus exp(i Jq Ldr) along all diffusion process trajectories, in our case along trajec- 
tories of Markovian process ([!]) with finite number of states. In this connection G(t) 
is the solution of equation 0: 

^- = PG- iLG. (3) 
dt v ' 

Here L is diagonal on stochastic variables Liouville's superoperator, which action 
on the spin transition operators B is determined by the rule: LB = H e B — BH g . 
Here H e,g = A e,g x (I e,9 ,m), I spin operator of a particle (for Fe 57 : I 9 = \,I e = 
|), rh - unit vector of particle magnetization direction, A e ' 9 - hyperfine interaction 
constants, P = P <8> E q - operator matrix diagonal on nuclear variables, E q - identity 
superoperator acting in quantum states space. 

Calculating integral (0) taking into account (|3|) one can get the following expres- 
sion for line shape @ 

N 

lp(lu) = ImJ2Pa{fFJ{a | L + iP- \E \~ l b)ffj j ), (4) 

o,6 
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where A = u + , V - natural half- width, fj - vector of polarization of falling 7- 
radiation, J - nuclear current operator, responsible for transitions between sublevels of 
ground \m g ) and excited \m e ) nucleus states, p a - relative probabilities of population 
of electron states a (there are N such states). So such calculation is reduced to 
spectral analysis of the generator of quantum evolution superoperator of a nucleus 
G~ 1 (X) = L + iP — XE , acting in space of operator-functions of the form: = 
\m e )(m g || a), which is a tensor product of the spin operators space and stochastic 
space of magnetic moment directions. Superoperator L, responsible for hyperfme 
decomposition of levels, acts by the rule 

I i I *') = [H e mem ,5 mgm , - H 9 m , mg 5 mem ,]5 aa ,. (5) 
Operator matrix P is diagonal on nuclei variables 

I P I = Paa'S me m' e Sm' g m g , (6) 

Paa = ~ ( fl I P I a ') ■ 

Relaxation matrix elements p aa > are transition probabilities in unit of time from state 
a to state a'. 

Of course one can consider the state |\& > at (|) and (^j) as usual state in direct 
product space of variables m e , m g , a and forget the superoperator nature of G — 
(L + iP - XE)- 1 . 

2 Spectral Analysis of Stochastic Matrices 

Let us suppose that easy magnetization axes coincide with axes of 3(4)-order of cube 
for the case of negative (positive) magnetic anisotropy constant. We accept the 
following indexing of these axes (fig. 1,2). For such indexing stochastic matrices of the 
Markovian process, giving magnetic moment motion, look in a following way 
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Spectral analysis of such matrices is well known [^|,P|,[10|. In particular P^ has two 
single eigenvalues and two triple ones (p = Ai = 0, r\ = A 2 = A 3 = A 4 = — 2(p+2g+r), 
£ = A 5 = A 6 = A 7 = — 4(p + q), ( = A 8 = — (6p + 4g)), and P^ has one single, one 
double and one triple eigenvalues (p = A x = 0, ( = A 2 = A 3 = — 6p, £ = A 4 = A 5 = 
Ag = — 4p — 2q). In connection with degeneration there is some arbitrariness in choice 
of eigenvalues, which can be eliminated by the requirement of definite symmetry 
of eigenvectors with respect to turns around quantization axis z (quantization axis 
direction for 3-axial case is different from || and this is connected with natural for 
giving case symmetry). In the case of N = 8(6) in the capacity of such eigenvectors 
let us choose states, that under turning at angle f (^f ) around quantization axis z 



obtain factor c n = i n ,n = 0, 1,2, 3(= exp(^f 1 
such vectors are the following 



n 



0, 1,2). In the case of N = 8 



K n = -[\1) + c n \2) + <£|3) + cl\A) + |1) + c n \2) + 4|3) + <£|4)], 
K n = \[\l) + c n \2) + 4|3) + 4|4) - |I) - c n \2) - 4|3) - 4|4)], 



and in the case N = 6 



[|l) + c n |2) + 4|3) + |I) + c n |2) + 4 



Kn 



|l) + c n |2) + C 2|3)-|l)-c n |2)-4| 



In this connection matrix P^ is diagonal in basis (Kq, K\, K 2 , K%, Kq, K\, K 2 , K3) 
and takes form 

p( 8 ) = ~diag(p,£,£,£,r),r],(,r)), 
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and matrix 

P^ = ^diag(p,C,C,U,Q 

in basic (K , K u K 2 , K Q , K h K 2 , ). 

It should be noted, that strictly speaking the use of probability transition matrices 
in produced form is certain assumption, i.e., is adjustment to real continuous diffu- 
sion process. However, in the case of small diffusion (or, that is the same, diffusion 
at strong fields) these matrices one can get directly from continuous diffusion model 
||,[|J on the base of strict asymptotic analysis of the diffusion operator. In addition 
it is found that nonzero transition probabilities are only those, that correspond to 
transition between potential wells, having attraction domains with common bound- 
ary. For matrix P^ this means that q = r = 0, and for P^ that q — 0. Moreover, 
calculations show that one can obtain symmetrical Markovian matrices only for the 
case of potential drift field. Non-symmetrical field automatically means that transi- 
tion probability from one easy magnetization direction to another is not equal to the 
inverse one. 

A 1 

3 Reduction of Superoperator G~ to Block Form 

Since superoperator G is invariant with respect to group C^{C^) action, it is nec- 
essary to study transformations properties of complete space of quantum-stochastic 
states respectively on such action, i.e., to construct basic elements of complete space 
of states, that under turning of coordinate system on angle f (^y) obtain factor c n . 
According to selection rules superoperator G in such basis is of block structure. 
For this aim it is naturally in quantum variables space to extract states with such 
properties relatively group C±{Cz) (such kind of states in stochastic variables space 
we have constructed above) and than to construct 4(3) nonoverlapping subspaces of 
complete quantum-stochastic space, the states of which obtain factor c n under turn- 
ing at angle f (y) around quantization axis z. How it is shown in ||, invariant with 
respect to group C4 subspace in the case N = 8 is divided into two nonoverlapping 
subspaces with different properties relatively time inversion (symmetry and asymme- 
try subspaces). Really, all others subspaces also suppose such division both in 4- and 
3-axial cases and also at other symmetries under proper disposition (consistent with 
the symmetry) of coordinate system. So it is possible to construct basis in which 
superoperator G consists of 8(6) equal blocks of order 8x8, that differ from each 
other only by diagonal elements, composed of stochastic matrix eigenvalues combina- 
tions. Nondiagonal part of each block is determined only by the transition structure 
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of nuclear system and does not depend on sign of magnetic anisotropy constant. 
Denoting by {n±| blocks ( s y mmetncal \ w jth respect to time reversal, that obtain fac- 



k asymmetric , 

tor c n = i n , n = 0,1,2,3 (exp(^p), n = 0,1,2) for 4(3)-axial case and denoting 
quantum transition operators 

^ = 1 2><2 U /?=I T )( T 1 ' 7=l 2 )( 2'' <T = 'l-><Tl' 

c ,lw-l, ,-lwl, ,3 —1 -3 1 

we have the following basic functions for blocks {n±} in 4-axial case 

\l) = K n (aT(3), \2)=K n+3 5±K n+1 e, |3) = K n+3l ± K n+1 a, 

|4) = K n+2 (n =F z/), |5) = i?„ +2 (/i ± i/), |6) = K n+3l =F ^n+i^, 
|7) = K n+3 5 t K n+1 e, |8) = K n (a ± 

and in 3-axial case 

\l) = K n (aTP), |2> = K n+2 5 ± K n+1 e, \3) = K n+2l ± K n+1 a, 

|4) = =F K n+2 u, |5) = i^n+i/U ± K n+2 u, |6) = i?„ +2 7 =F K n+l a, 

|7) = iT n+2 5 =F K n+1 e, |8) = ± 0), 

Here in the 4(3)-axial cases all indexes are taken by modA(3). In this basis diagonal 
submatrix P n ± takes form for 4-axial case 

P 0± = -diag(A, n, n, n, E, T, T, T) (7) 

. n n n 

Pi,s± = -diag(Q, —, —, — , T, 6, 9, T) 

P 2 ± = —diag(Q, Q, fl, A, T, T, T, S), 

where A = 0, Q = 4p + 4q, E = Qp + 2r, T = 2p + 4q + 2r, 9 = 4p + 2q + 2r. For 
iV = 6we have 

p 0± = -diag(A, n, n, n, r, r, r, r), (8) 
, n n n 

Pi,2± = -diag(Q, —, —, — , T, T, T, T), 
where A = 0, f2 = 6p, T = Ap + 2q. 
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Superoperator L action at all received subspaces independently on symmetry ap- 
pears to be the same and is given by matrix 



L 



( D 
ll>t o 



D 



J 



A e /y/2 (2A e - A 9 )/y/6 (A e - A 9 )/2V3 \ 

A e /y/2 (A e + A9)/2y/3 -(A 9 + 2A e )/y/6 

-Av/VQ (3A e -A 9 )/2V3 A e /V2 

\(3A e + An)/2y/3 -A 9 /VE A e /V2 

Block of form D was introduced first in for invariant subspace of quantum- 
stochastic states antisymmetrical with respect to time inversion (we denote this sub- 
space as 0-) in the case N = 8. Remarkably, that other blocks of superoperator L 
(that have another symmetry, i.e., that obtain factor c n under turning of coordinate 
system on angle | (^f)) have the same form. Moreover, one can show that the same 
form of elementary block remains for certain other symmetries, for example for tetra- 
hedron symmetry, cylindrical symmetry and some others. So elementary block D is 
in some sense universal. 



4 Spectral Analysis of Quantum-Stochastic Super- 
operator 

Last paragraph is devoted to reduction of quantum-stochastic superoperator to block- 
diagonal form convenient for calculations its spectral characteristics (eigenvalues and 
eigenvectors). This problem can be solved under assumption that one of terms L or 
P composite superoperator G" 1 is small in comparison with another. It is naturally 
to assume smallness of stochastic term P. Note that eigenvalues of unperturbed 
superoperator L are known explicitly. These eigenvalues are the super-fine transition 
frequencies X ±i = ±b h where &i = \{3A e - A 9 ) , b 2 = \{A e - A 9 ) , 6 3 = -\(A e + A 9 ) , 
64 = — |(3v4 e + A 9 ) . Corresponding eigenfunctions have the form [[| 

* ±1 = + + -L|3) + ^|4) ± (i|5) + f |6) + -L|7> + ^|8» 

** = > - 2> - >> + ± <">> ~ 5> + l m + 0) 
= i> + > + - > ] ± ( 5> - > + ° + > 



s 



- > } + ^> - > ± ( x' 5 > 4' 6 > + ^ - >»■ 

Standard perturbation theory in the case of absence of degeneration let us calculate 
eigenvalues corrections 

AA fc = (tf fc | iP | tf fc ) 

and also eigenvectors corrections. In our case for eigenvalues this leads to results 
reduced in Tab.l of eigenvalues corrections. 



Tab, la. Eigenvalues corrections for 4-axis case. 





block 0± 


block 1,3± 


block 2± 


6X1 


i(19/9p + 2Q/9q + r) 


i(31/9p + 26/9g + r) 


i{3p+ 10/3g + r) 


5X 2 


i(3p + 8/3q + r) 


i(3p + 8/3q + r) 


i(3p + 4q + r) 


sx 3 


i(3p+10/3q + r) 


i(3p+10/3q + r) 


i(3p + 2g + r) 


5A 4 


i(35/9p + 28/9q + r) 


i(23/9p + 28/9q + r) 


i(3p + 8/3g + r) 


Tab.] 


b. Eigenvalues corrections for 3- axis case. 





block 0± 


block 1, 2± 




i(3, 5p + q) 


i(4,25p + g) 


5A 2 


i(Ap + g) 


i(4p + 9) 


5A 3 


i(4, 5p + g) 


i(3, 75p + q) 


5A 4 


i(5p + g) 


i(3,5p + q) 



5 J-Selfadjointness and Qualitative Spectrum Be- 
havior 

Reductions of operator matrix L + iP to block-diagonal form at mentioned basis 
permit to receive effective way of calculating of line shape reduce to inversion of 
matrices of 8 x 8-order corresponding to each block. In the case of equilibrium initial 
distribution resonance absorption spectrum is described by the only block of the 
superoperator L+iP, namely by that, which is invariant with respect to turns on angle 
§(^-) and antisymmetrical relatively time inversion. But if the initial distribution is 
arbitrary one it is necessary to know all blocks of the superoperator L + iP for line 
shape determination. 

The superoperator L + iP spectrum is, generally speaking, complex one, never- 
theless, it is possible to get qualitative information on spectrum behavior at different 
relations between parameters of magnetic and electron spin systems. The superoper- 
ator L + iP spectrum lies at lower half-plane and is placed symmetrically relatively 
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imaginary axis. This follows from J-selfadjointness of L + iP. Really 



when J=idiag(l, 1, 1, 1, —1, —1, —1, —1). 

Let us consider for illustration the case of equilibrium initial distribution, which 
is described by block 0-. At the limit case of fast relaxation (at high temperatures) 
eigenvalues of this block (just as of any other block) belong to imaginary axis, and 
this case corresponds to resonance absorption spectrum consisting of the only line of 
Lorentz width on the frequency to. At diminishing of relaxation eigenvalues as before 
some time belong to imaginary axis moving on it, because at fast relaxation limit there 
are two single and two triple eigenvalues both in 4- and 3-axial case (0),©, but due to 
the J-selfadjointness eigenvalues can leave imaginary axis only by pairs. At certain 
moment one of pairs will diverge (and also this can be not one pair from initially 
existing pairs in triple eigenvalues but one of anew formed in a process of random 
walk along imaginary axis). In this connection resonance absorption spectrum will 
have a form of two lines of certain width, symmetrical relatively frequency uj. Under 
further diminishing of relaxation parameters new pairs of lines appear, and the width 
is decreased. At limit all eigenvalues of the matrix come down from imaginary axis 
and are dropped on real axis, that on resonance absorption spectrum correspond to 
line shape with peaks on frequencies uj + b±i (frequencies 6±4 are prohibited and don't 
appear on spectrum), the width is decreased as far as eigenvalues approach real axis, 
i.e., under diminishing of relaxation or, that is the same, of temperature. 

It is necessary to note that symmetrical line shape corresponds to potential drift 
field (in DOM — to symmetrical matrices P). Calculations using non-symmetrical 
transition probabilities demonstrate non-symmetrical line shape with respect to main 
frequency uj. 

Author is grateful to K.A.Makarov and B.S.Pavlov for useful discussions. This 
work is supported by RFBR grant 00-01-00480. 
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fig. 1. Mutual disposition of the easy magnetization directions and 
coordinate system (4-axis case, N=8). 

Values of unit vector rh of particle magnetization 
rh = ±rhi, % = 1,2, 3, 4. 

mi = (y/2/3, 0, 1/V3), m 2 = (0, y/2/3, 1/^3), 
m 3 = (-y^A 0, 1/V3), m 4 = (0, -y^tt, 1/^3). 

^asy magnetization directions coincide with ones from the center of cube to its 
vertices. 
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fig. 2. Mutual disposition of the easy magnetization directions and 
coordinate system (3-axis case, N=6). 

Values of unit vector rh of particle magnetization 
rh = ±rhi, i — 1, 2, 3. 

Eas?/ magnetization directions coincide with ones from the center of cube to the 
centers of its faces. 
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